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The fluctuations of capacitance of a two-dimensional island are studied in the regime of low 
electron concentration and strong disorder, when electrons can be considered classical particles. The 
universal capacitance distribution is found, with the dispersion being of the order of the average. 
This distribution is shown to be closely related to the shape of the Coulomb gap in the one-electron 
density of states of the island. Behavior of the the capacitance fluctuations near the metal - insulator 
transition is discussed. 

PACS numbers; 71.55.Jv 



Capacitance C is conventionally understood as a well 
defined geometrical property of a metallic sample. For 
example, for a metallic sphere of radius R capacitance 
C = R, and when the sphere is large, in the first approx- 
imation capacitance does not depend on the distribution 
of impurities inside the sphere or on its charge. How- 
ever, in very small metallic samples fluctuations of ca- 
pacitance become observable. Recently, such fluctuations 
were measured in the semiconductor quantum dots as a 
function of the total charge of the dot using the Coulomb 
blockade phenomenonHn. In such an experiment a small 
quantum dot is weakly coupled to current leads while a 
gate is placed in the proximity of the dot and is used 
to vary its electrostatic potential. At low temperatures 
the charge of the dot is typically quantized and there is 
no significant current. However the gate voltage can be 
tuned in such a way that the ground states with N and 
A'' -I- 1 electrons are degenerate. At this gate voltage cur- 
rent can flow through the dot. The resulting conductance 
vs. gate voltage comprises a series of sharp peaks (charg- 
ing spectrum). The spacing between two peaks AVg can 
be expressed in terms of the ground state energies Ef^ of 
the dot with A^ electrons: 
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Here a is the geometrical coefficient, Cn is the capac- 
itance of the dot with A^ electrons. This equation may 
be considered as a definition of the capacitance. For a 
macroscopic body with the positive background charge 
eA^Oj En has a simple form: En = e'^{N — NqY /2C and 
Eq. (|l|) gives Cm = C = const. For the quantum dot 
the charging energy A jv was foundcl to have surprisingly 
large relative fluctuations: 
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Here < ... > denotes the averaging over A^. Much ef- 
fort has been done to explain such large fluctuations. 
First, the experimental data were compared with the 
so-called constant interaction model in which Ajv = 
e^/C + r/pf+i — rjN, where r/pf is the A^-th single-electron 



energyH. The fluctuations of the spacing between the 
nearest neighbor levels are of the order of the average 
spacing E-p /N, where Ep is the Fermi energy. Hence, for 
a dot of radius R 
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where = as is the screening radius of the two- 
dimensional electron gas, qb is a semiconductor Bohr ra- 
dius, which is close to lOnm in GaAs. For this value of gb 
and for R - 200nm (see Ref. |[) Eq. (|) gives fluctuations 
that are substantially smaller than observed in the ex- 
periment. This discrepancy initiated the computer mod- 
elling and the analytical calculations of 6 (see Ref ||,||). 
When discussing theoretical results one should keep in 
mind that in all the theoretical works S was obtained 
by averaging over the different realizations of disorder, 
instead of number of electrons A^. Below we will also 
assume that in strongly disordered systems there is no 
difference between these two definitions of fluctuations. 

Analytical diagrammatic calculations based on RPA 
confirm Eq. (p[)t^_pn the other hand, the results of com- 
puter modellingElB agree with Eq. (jsj) for weak interac- 
tions (large Ts and as) and lead to larger and interaction 
independent S for strong interactions, corresponding to 
low density electron gas. This fact was identifiedpas a 
failure of RPA in the the low density electron gasQ. It 
can be interpreted easily in terms of revision of equation 
''s = Qs at low densities n <^ a^^. Indeed cannot 
be smaller than the average distance between electrons 
n~^/^ and at the small densities one should substitute 
Ts — n~^/^ into Eq. (||). It is not clear yet whether 
such a simple modification of Eq. (^ caujfluantitatively 
explain numerical and experimental datao, which seem 
to indicate that S is almost R independent. Hence, it 
is challenging to understand what happens with 6(R) in 
the limit of a very low electron density. 

In this paper we theoretically study the fluctuations of 
capacitance of the island in the extreme classical limit 
when the quantum kinetic energy of electrons is much 
smaller than both the disorder strength and Coulomb in- 
teractions. We consider the case of a large disorder when 
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the ground state of the island is a Coulomb glassQ. Below 
we show that for a piece of Coulomb glass or, in other 
words, a Coulomb glass island, the fluctuations are large, 
S is of the order of unity, does not depend on R, i.e. is 
universal for a given shape of the island. For the square 
sample we find S = 0.32. (Previously a similar statement 
about giant, of the order of unity, relative fluctuations of 
the polarizability of the Coulomb glass island was made 
in Ref. |). 

We study probability density of Ajv and show that it is 
a universal function of the ratio x = A^/ < Aat >. We 
also discuss how the transition from the Coulomb glass 
to metal is reflected in the capacitance fluctuations. 

When an electron is added to a large metallic sam- 
ple its charge is distributed in the unique way according 
to the electrostatic theory. For this reason with addi- 
tion of every electron the electric potential grows by the 
same amount e^/C, or in other words C — const. In 
the Coulomb glass the electronic states are localized, so 
that every new electron is put into some localized site. 
Then electrons rearrange themselves in the vicinity of 
this site. However it was shown in Ref. ^ that such rear- 
rangement happens only with probability close to 1/2 on 
every scale. As a result the added charge in the majority 
of cases is localized in the region smaller than the size 
of the island. When next electron enters the island its 
charge is centered near another site in the island. The 
distance between this site and the position of the previ- 
ous electron fluctuates between and 2R. As a result 
the difference between energies required to bring two se- 
quential electrons. An and capacitance Cn experience 
roughly speaking hundred percent fluctuations. 

To verify this reasoning we study the capacitance fluc- 
tuations numerically. We use thcrjlattice model of the 
Coulomb glass suggested by EfrosQ and widely used to 
study the Coulomb gap in the density of states (DOS). 
The Coulomb glass island is modelled by the square 
M X M lattice, with every site being either empty ( oc- 
cupation number = ) or occupied by one electron 
(rii = 1). Electrons interact with each other by Coulomb 
interaction. The interaction energy between the near- 
est lattice sites is chosen to be the unit of energy and 
the lattice constant is the unit of length. Disorder is 
introduced by the random site energies which are dis- 
tributed uniformly between —1 and 1. The corresponding 
Hamiltonian has the form 
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Here Ui is the potential due to the uniform background 
charge making the system electrically neutral for N elec- 
trons. We find the ground state energies of TV— 1, N, N+1 
electrons where N is the integer part of M^/2 and then 
calculate A at using Eq. (^. To find the ground state 
we use two different methods: the exhaustive enumera- 
tion and the simulated annealing. The first one is used 
for relatively small samples with M < 5. We enumerate 



all the possible states of N electrons on M x M lattice 
sites and find one with the lowest energy. In the second 
method we employ the simulated annealing technique, 
running the finite-temperature classical Monte-Carlo for 
some time and taking the lowest energy state. The con- 
vergence of the solution to the ground state has been 
checked by doubling the time of the simulation and mak- 
ing sure that the solution is not affected. The reliable 
results have been obtained by this method for M < 8. 
Atv has been calculated typically for 1000-2000 different 
realizations of disorder. We have obtained 5 ~ 0.32. Nor- 
malized distributions F{x) of the ratio a; = Ajv/<Ajv> 
for M — 4,5,7,8 are shown in Fig. [I]. We have found 
< Ajv >= 2.3/ M, in a good agreement with the inverse 
capacitance of metallic square of the same size. Remark- 
ably, F{x) does not depend on M. We emphasize also 
that contrary to the predictions of the constant interac- 
tion model the inverse capacitance of the Coulomb glass 
island can be both larger and smaller than that of the 
metallic island of the same shape. 
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FIG. 1. The inverse capacitance distribution is pre- 
sented for different sample sizes. The line is the fit by Eq. 

(i- 

The fit in Fig. |l| is given by the following equation: 

= ( ^ / o/ ^2^ (5) 

^ ^ \ 4(j;-xo)exp(-2(a:-a:o)^) , x>Xo, ^ ^ 

where xq = 0.37. Two interesting features of F (x) are 
clearly seen. Firstly, this function has a termination 
point at X = xq. One can easily check that Ajv cor- 
responding to this point is equal to the smallest possible 
Coulomb interaction : l/r,„ax, where Tmax = {M — 1) 
is the maximum distance between sites in the square 
M X M. Secondly, F{x + xq) is identical to the Wigner 
surmise for the nearest neighbor distance distribution of 
the levels of a random matrix. We shall show below that 
this is only an interesting coincidence. 

Now we will interpret both features establishing the 
relation between F{x) and the one-electron DOS g{e) of 
the Coulomb glass island. The energy of the one-electron 
excitation localized at i-th site can be written as 
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For an empty site, for example, it is the energy required 
to bring an electron from infinity to this site. When aver- 
aging this DOS over different realizations of the disorder 
potential or, in other words, over different samples 
match the chemical potential in them {fj, - averaging)! 
The chemical potential of the island is situated halfway 
between the largest energy of the occupied states and 
the lowest energy of the empty ones. The corresponding 
DOS is shown in Figure ^ 




and capacitance. The excitations relevant tn-jcapacitance 
are electronic polarons introduced by EfrosQ. Their en- 
ergies ii can be used to calculate A^. Indeed, ii is de- 
fined as the energy required to bring an electron to the 
site i and rearrange the other electrons in order to reach 
minimum of the total energy. From this definition it is 
obvious that the minimum polaron energy of the empty 
states is i?Ar+i — i?Ar, and the maximum polaron energy 
of the occupied states is En — E^-i- It is clear that Ajv 
is just the energy gap between these two polaron states. 
It is known, however, that in two dimensions thCjBolaron 
energies are very close to the one-electron oneaa. This 
means that A^v can be well approximated by the energy 
difference between the lowest empty and the highest oc- 
cupied one - electron states. This immediately explains 
the existence of the termination point of F{x): A^r can- 
not be smaller than the smallest interaction between two 
electrons within the island. Moreover the function F{x) 
can be related to the one-electron DOS g{e) in the fol- 
lowing way. 

As the inverse capacitance of an island is equal to 
the difference between the lowest empty and the high- 
est occupied states' energies, our problem is to find the 
distribution function of this difference. The probability 
P (A > e) to have it bigger or equal than certain value e 
is equal to the probability not to find energy levels in the 
region —I < e' < |. Assuming the Poissonian statistics 
of the level distribution we arrive at: 



FIG. 2. DOS for the 8x8 sample averaged over disor- 
der. The main picture shows the region near the Fermi 
level. The Coulomb gap in the density of states for an in- 
finitely large sample [Eq. (Q)] is presented by the straight 
lines. The inset shows the general view of the DOS. 

The important feature of this DOS is a linear Coulomb 
gap, which at the small energies crosses over to the hard 
gap related to the finite size effects □. Linear dependence 
of the DOS for |e| > 0.2 agrees with the analytical ex- 
pression for an infinitely large sample 
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derived in Ref. ^,0. The total width of the hard gap is 
equal to 1/rmaxi where r^ax is the maximum distance be- 
tween two points in the island. Indeed, the energy that 
is required to transfer an electron from site i to an empty 
site j is equal to: 
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The minimum difference between the energies of empty 
and occupied states cannot exceed the minimum interac- 
tion energy within the island and, hence, is greater than 
or equal to l/r^ax- 

Let us now explain how the one-electron DOS can be 
used to find F{x). Strictly speaking the one-electron en- 
ergies are not directly related to the ground state energies 
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Here g(e) is assumed to be normalized to unity: 
g{e')de' = 1. The expression in the exponential is 
the average number of electrons in the band of energies 
(— |, I). The probabihty density of A is, hence, equal to 
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This expression establishes the general relationship be- 
tween the one-electron DOS and the inverse capacitance 
distribution function for the Coulomb glass island. Hav- 
ing been applied to the numerically obtained DOS it gives 
a very good agreement with the actual F{e) . Substituting 
Eq. @) into Eq. one arrives at the Gaussian asymp- 
totic behavior of F{x): lnF{x) cx —x^, 1 ^ a; <C M. 
This asymptotic, however, does not persist in three di- 
mensions, where the DOS has a different form. 

Let us now examine what happens with d when a hop- 
ping term is added to the classical Hamiltonian, given by 
Eq. i): 
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Here J is the hopping matrix element and a] is the cre- 
ation operator of the electron at site i. The summation is 
carried over the neighboring sites. The ground states of 
this Hamiltonian were found numerically using the Lanc- 
zos algorithm. The system considered was of size 4x4 
lattice sites, with up to 7, 8 and 9 electrons. The results 
are depicted in Fig. ||. 
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FIG. 3. The relative inverse capacitance fluctuations 
of 4 X 4 island as a function of the hopping matrix element 
J. 

As it is seen that 6 decreases by a factor of 2 from J = 
to J = 0.4. Such a decay is consistent with the tendency 
of metallic samples to have smaller capacitance fluctua- 
tions. At the same time the shape of F{x) becomes more 
symmetric in agreement with Ref. Unfortunately we 
were not able to do such calculations for M > 4. Hence, 
the existing numerical data leave open the challenging 
question of how the crossover happens between S ^ 1 
in the classical case and Eq. in the quantum one for 
large enough samples. Below we try to answer this ques- 
tion concentrating on the three-dimensional case where 
the insulator-metal transition happens at some critical 
value J — Jc- We assume that this transition is accom- 
panied by the divergency of the wavefunction correlation 
length ^ and the dielectric constant k: ^ = \Jc — J\^'^ 
and K = I Jc — J|~^- It was argued that k ^ (^/r^)^, 
where fs is the screening radius of the three-dimensional 
degenerate Fermi-gas, or C = 2i/ ( see Ref. || ). When 
J approaches Jc from the insulator side the growth of k 
plays a very important role in the distribution of charge 
of the added electron, even when ^ is still much smaller 
than a sample size R. It is well known that if a localized 
charge is put inside a dielectric sample with k ^ 1 the 
sample becomes polarized in such a way that almost all 
of the added charge appears on its surface in the form 
of induced charge. Only a small fraction of charge c/k 
remains localized inside. If the second electron is added 
to the sample at the distance r from the first the in- 
teraction energy between these two consecutively added 
electrons fluctuates by b^/kR ^ e^/R, provided that r 
fluctuates in the range < r < 2R for a sphere of radius 
R. Therefore 

(r,/C)' - (Jc - Jf" (12) 



Eq. (|T|) is valid only if ^ < i?, or J < Jc - A J, where 
AJ = {rs/R)^^^ ■ At = i? the relative capacitance 
fluctuations saturate at 

5--{rslRf (13) 

Eq. ( p^ ) is the three-dimensional analog of Eq. (||). It can 
also be obtained from the assumption that fluctuations of 
Ajv are equal to the fluctuations of the spacing between 
the one-electron quantum levels. Predicted behavior of 
function 5[J) is schematically depicted in Fig. ^. 
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FIG. 4. The relative inverse capacitance fluctuations 
of the three-dimensional island as given by Eq. ( |l^ ) for 
J <Jc~AJ and by Eq (|l|) for J > Jc - AJ. 

We conjecture that in the two-dimensional case, where 
as it is commonly believed ^ grows monotonically with 
J, the crossover from (5 ~ 1 to Eq. (||) happens in the 
similar way: S ^ r^/^ at ^ ^ ^ i? and S ^ rg/R at 
^ R. Schematically this behavior is similar to the one 
shown in Fig. ^. Our numerical results do not contradict 
to this prediction. 

Up to now we have dealt with the unscreened Coulomb 
interaction between electrons. In case if a metallic gate 
is situated at a small distance d from the plane of the 
island the long-range part of the Coulomb interaction 
is screened. Let us now discuss how this screening af- 
fects our results. Consider first the case of an extremely 
close gate when one can completely ignore the interac- 
tions between electrons. In this case the constant (zero) 
interaction model is valid, N electrons occupy N lowest 
one-electron levels ry^, and Ajv = Vn+i—Vn is just a near- 
est neighbor level spacing. At J = it has the Poisson 
distribution, therefore i5 = 1. In the metallic phase (large 
J) the random matrix theory is applicable and S = 0.52 
(see Ref. ^). We studied this crossover for up to 10^ real- 
izations of disorder in the square 4x4, with the Coulomb 
interaction replaced by 1/r - l/{r^ + {2dfy^^, d = 1/2. 
We obtained S = 0.67 at J = and (5 = 0.39 at J = 0.4 
in the qualitative agreement with our expectations. At 
J = we also observed a drastic change in the form of 
F{x): the gap at small positive x disappears and F{x) ap- 
proaches the Poisson distribution. The most remarkable 
feature of the observed distribution is that there exists a 
very small tail of F{x) at x < 0, so that the differential 
capacitance of a strongly screened disordered island can 
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be negative. The possibility of the negative capacitance 
in the presence of a gate was recently demonstrated in 
Ref. H. The observed probability to have negative capaci- 
tance is consistent with the calculations of these authors. 
For our set of parameters it is equal to 3 • 10^^. 

In conclusion, we have found the large universal rela- 
tive fluctuations of capacitance of the Coulomb glass is- 
land and described a scenario of their decay when the 
system undergoes the insulator-metal transition. We 
are grateful to M. M. Fogler for valuable discussions. 
This work was supported by the NSF Grant No. DMR- 
9321417 and by University of Minnesota Supercomputer 
Institute. 
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